Isolated and dynamical horizons from a common perspective 

Mikolaj Korzynski 
Institute of Theoretical Physics, Warsaw Universit^ 
(Dated: February 7, 2008) 

Abstract 

A framework is developed in which one can write down the constraint equations on a three- 

I dimensional hypersurface of arbitrary signature. It is then applied to isolated and dynamical 

O ' 

' horizons. The derived equations can be used to extract physicaly relevant quantities describing 
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the horizon irrespective to whether it is isolated (null) or dynamical at a given instant of time. 
Furthermore, small perturbation of isolated horizons are considered, and finally a family of axially- 
symmetric exact solution of the constraint equations on a dynamical horizon is presented. 
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I. INTRODUCTION 



Since their introduction, isolated and dynamical horizons have been used to describe 
Mac. Kole. i„ gene.a, dynamical situation. both analytic |1 B fl fl B and numerical 



context. Both are defined as three-dimensional hypersurfaces foliated by mariginally 
trapped surfaces (so called mariginally trapped tubes, MTT's). They are supposed to repre- 
sent the boundary of a black hole in a non-stationary spacetime. Isolated horizons (IH) are 
null hypersurfaces and describe black holes which are in equilibrium with the environment 
at the moment. Dynamical horizons (DH) on the other hand are spacelike and describe a 
black hole absorbing matter or gravitational radiation. 

In physical situations we expect an isolated horizon to become dynamical if matter or 
gravitational radiation falls on it. On the other hand, it is believed that a perturbed black 
hole will return to a stationary state after the matter infiow ends and all gravitational 
radiation is radiated away. In this case a DH should asymptotically become null. In both 
processes we are led to investigate a hypersurface which can change its signature (from null 
to spacelike) depending on circumstances. 

However, geometrically null and spacelike surfaces are very different. The latter are 
equipped with a well-defined, positive definite metric, which makes them Riemannian man- 
ifolds. Such hypersurfaces are completely described by their metrics and extrinsic curvature. 
On the other hand, on null surfaces one usually uses a different set of geometrical objects 

fl. 

It is tempting to use perturbative techniques to investigate small perturbations of an 
IH IllL Il2l | . This is not easy, as in presence of smallest perturbations the IH is not 

strictly null anymore. On the other hand, the ortogonal vector cannot be normalized if a 
DH becomes null and therefore the standard definitions of extrinsic curvature and covariant 
deriavative do not make sense any more j^. 

This demonstrates the need to introduce a way of description of a MTT (or any hyper- 
surface) which is applicable to any signature. It should involve "signature-blind" variables 
which do not refer to any normalization of its ortogonal vector or to the three-dimensional 
metric on the horizon. It should however contain all information needed to recover the 
standard geometrical objects if the surface happens to be spacelike. 



Although variables of that kind have been used in calculations by several authors ll 
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there hasn't been a complete discussion of the Einstein constraint equations in this 
context, neither on MTT's, nor in general. The aim of this paper is to introduce such 
variables and discuss the Einstein constraint equations on any hypersurface - spacelike, null 
or timelike. We then focus on MTT's. We also analyse a lineary perturbed IH and finally 
present a family of axialy symmeric, exact solutions of the Einstein constraint equations. 



II. "SIGNATURE BLIND" VARIABLES 



We consider a 3-dimensional hypersurface H foliated by a family of 2-dimensional space- 
like surfaces of spherical topology, called sections. The induced metric on H may have 
arbitrary signature, or even be degenerate. 

By convention, the greek indeces will denote spacetime objects, while the latin indeces 
objects defined on sections. In particular, if is a one-form defined on spacetime, 
will denote its puUback to a section, whereas if X° is a spacetime vector, X"^ will be the 
ortogonal projection to the tangent space to sections. 

On each section we have a positive definite metric qab- At each point we can pick up two 
null vectors, fc^ and l^, ortogonal to the leaves of the foliation and normalized fc^ = — 1. 
We also choose a vector field tangent to H, preserving the foliation and ortogonal to 
sections. 

We may now decompose n according to = nil'^ + fc^. The vector field = 
rii — rik is ortogonal to both leaves of the foliation and n, and therefore ortogonal to H. 

The choice of I and k is not unique, as we may rescale them simultanously by any non- 
vanishing function A: l'^ = Al'^, k'^ = A^^ k^. This gauge freedom can be partially fixed 
by first fixing / and k at one section and then imposing the conditions 

V„/^ = -V^l,k>^ 

Vnk'^L = Vnk''k^ = 0. (1) 

They can be checked to preserve the normalization of / and k, their null character and 
their ortogonality to the foliation (see Appendix [X]) . 
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We have therefore specified the n component of the covariant derivative of / and k. Its 
pullback to sections constitutes so called Weingarten tensors or deformation tensors 

VaIb = Lab = o^'ab) + (^) 
VAks = Kab = cr^AB) + Ub- (3) 

They have been decomposed to traceless, symmetric part, called shear, and the trace called 
expansion. The antisymmetric part vanished since k and I are surface-forming. 
We also have the rotation one-form uoa = — Va^'^ k^ = V Ak^ ^ti- 

If H is spacelike, the data specified here are equivalent to the standard objects in con- 
straint equations, i.e. the three-dimensional metric and the extrinsic curvature. The precise 
transformation is given in Appendix ^ 

We still have the gauge freedom to multiply and k^ by one function at a given section: 

/'M = FF 
k'^' = p-^k" 

uj'a = uJA-DAln\F\ (4) 



III. THE EINSTEIN CONSTRAINT EQUATIONS 

The constraint equations on any three-dimensional hypersurface are the Einstein equa- 
tions involving the Gq/jt" component of the Einstein tensor, where r is ortogonal to the 
surface. 

The Gq/jt" one-form can be decomposed to the pullback to sections GaB and two 
sections-ortogonal components Ga/3T°'l'^ and GapT^k^. 

We shall denote the covartiant derivative on the section 2-spheres, compatible with g^s, 
by D, its Laplacian q"^^ DaDb by A and its Ricci scalar by 71. The constraint equations 
read 

Cn 00 A = G^AT^ - ni{DB L\ + 00 A O') - DbUi L\ + 

+nkiDBK'' A-iOAe'') + DBnkK''A (5) 
Cji = -G^^T^F + Auk - 2cj^ DAUk -nk(^ + Da uj^ - uauj^ + 9^ 9''^ + 

-uiLabL^'' (6) 
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G.^T^'k" + Am + 2a;^ DaUi - m { ^ - Dau;"" - uaoo'' + 6'' ] + 



-Uk Kab K 



AB 



(7) 



We have introduced the notation Kab = o-ab ~ Qab and analogically for Lab- These 
equations, together with the evolution equation for the metric qab 



Qab = 2ni Lab + 2nk Kab 



(8) 



constitute a system, which yields initial data for Cauchy problem in the spacelike case. It 
can be analysed and solved, as we shall see, at least in some special cases. 

If we introduce a "time" coordinate on H compatible with n and the foliation, i.e. t is 
constant on sections and n'^t,^^ = 1, we can replace the Lie derivatives in ((SI)-® by ordinary 
partial derivatives with respect to t. 



Equation (0) was discussed in 
Equations (jH)) and ((Tj) appeared in 



81 in the null case and 
151 and 



0. 



in Q i 



in a general situation. 



IV. CONSTRAINT EQUATIONS ON A MTT 

We now set 9i = 0. The equations ©-((HI) become now 

^nQAB = '^ni (Jab + 2/2^ Kab (9) 

Cn UJA = G^AT^ - DBini a' \) + nu{DB K^a-uja + DbUu K^'a (10) 

= -G^.r^/" + Pnu-ni a^^ (11) 

CnQ"" = G.^T'^k'' + P*ni- nu [a\B + \o' ') • (12) 

We have introduced two elliptic differential operators 

P = A - 2uj^ Da - + div iu - cu^j (13) 

P* = A + 2cu^ Da-(^- div uo-uo^Y (14) 



They are adjoint to each other [15|: for any section S and complex functions f,gonS we 
have 

/ fPg^/deUid^x= j g P* f ^/detq d^x. (15) 

'J S " s 
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Given the matter sources (represented by the Einstein tensor), we have a system of 7 
differential equations for 12 quantities. All but one are evolution equations. Equations (fTT|) 
and (jl2|) appeared in 15 1 and their combination was used in 



Equation (jlip doesn't involve time derivatives. It is a differential equation for and 
algebraic for rii. Given a spacetime and a mariginally trapped surface therein, we can use 
pi|) to construct a family of mariginally trapped surfaces constituting a horizon. We may 
solve it by specifying Uk and then deriving rii or, if P is invertible, the other way round. In 
both cases we have a freedom of specifying one function. 

The character of H is determined by t'^ = 2nink- It is null if = or rii = 0, spacelike 
if the signs of rii and are opposite and timelike if the signs are the same. 

If cr' = and G^^i, = 0, we can take = and any rii. In this way we obtain a 
non-expanding horizon. Although we have a complete freedom of specifying one, positive 
function n^, various choices do not lead to different submanifolds of the spacetime, but 
rather change the slicing of a given H. With proper choice of slicing we can make H a 
weakly isolated horizon. 

If any of those quantities does not vanish and P turns out to be strictly negative, one can 
prove that H is strictly spacelike provided that rii > 0. This can be done by applying the 
strong maximum principle |[16|]. In simple terms, one proves that P~^f is negative if / > 0. 
Note that although P itself is not invariant with respect to gauge transformations (@]), the 
property stated above holds even after rescaling of I and fc, provided that we choose F > 
(we don't flip the direction the vectors point). 

In spacelike case different solutions of (fTT|l lead to truly different horizons, i.e. different 
submanifolds of the spacetime [l^]- Furthermore, if 9^ < 0, the hypersurface is a DH. In 
17[ it is also proved that, under some reasonable technical assumptions, there are strong 
restriction concerining the location of DH's with respect to each other. In particular, they 
cannot foliate any neighbourhood of DH. 

A. Angular momentum 

Following Q and 0, 

we define the angular momentum of H with respect to a specified 
vector field X"^, tangent to a section, as 

Jx = / X^iUA Jd^qd^x. (16) 

ottG Js 
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In the null case this definition agrees completely with the standard one. In the DH case it 
is slightly different from the definition in [l|, which in our terms can be written as 



J'x 



1 



SttG Js 



uja + --Da In 



det q (fx. 



(17) 



However, as was pointed out in jlj , the definition is not reliable if is not a symmetry, or 
at least divergence-free. If DaX^ = 0, both definitions actually coincide. Moreover, (fT^ is 
gauge invariant with respect to (jH). 



B. Small perturbations of an IH 

The framework introduced above is particularly well-suited to discuss the perturbation 
of an IH. As we mentioned, an IH can be perturbed by matter infiow or (T^^r" ^ 0) or by 
external gravitational field {cxi 7^ 0). In this paper we will analyse the perturbation by a 
gravitational field in the linear order. 

External gravitational field, for example a gravitational wave, causes a' 7^ 0. This affects 
the equations (fTUj) and ©, but not (fTT|) . as it involves cr' ^. Therefore in the linear order Uk 
remains 0. 

This has important consequences: first, the second term in (jH)) vanishes and therefore 
the metric is perturbed only by riiai. This means that even if the geometry of sections is 
deformed (by tidal forces), the deformation does not affect their total area. 

Moreover, one can integrate equations and (|Tn|) with the Uk terms absent. They yield 

qAB{t) = qAB{to) + f 2ni{t') a\j,{t') dt' + o{a') (18) 
J to 

uJA{t) = coAito) + D^'iqABit) - qAsito)) + o((t') (19) 

[D^ in the second equation is take with respect to qAB{to))- This means that while u (and 
therefore the angular momentum of the horizon) may change in the result of influence of 
external gravitational fleld, its change is rigidly connected with the change of geometry. 
If the horizon returns to its original shape after the perturbation, the angular momentum 
returns to its original value. 

The conclusions above are valid for, among other things. Booth and Fairhurst's slowly 
evolving horizons 0, Q| without matter accretion. Under their assuptions, cr' is of the order 
of eajj Is], e being the Booth-Fairhurst "small parameter" and an denoting the horizon 
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area. However, we do not make any use of the slow change of 9^ 



(condition 3 in their papers). 

V. A FAMILY OF EXACT SOLUTIONS 

Bartnik and Isenberg have investigated the constraint equations on a DH in the case 
of full spherical symmetry We can generalize their result to a special class of axialy 
symmetric horizons. 

The set of equations (j^- (fTT|) seems very difficult to treat at first glance. However, we 
may achieve a great deal of simplification if we assume rii, n^, 9i, 9k, T^^t'^I^ and T^pr'^k^ 
to be constant on sections. In this manner we hope to reduce the equations (P)- (fTT|) to a 
system of ODE's. For simplicity we also put = (Ji = 0. This implies, due to 0, that the 
horizon geometry simply scales by a constant conformal factor with time q{t) = e^^*-* go- 

The assumption of scalar quantities to be constant on sections imposes further conditions 
on the geometry of H. In order to satisfy both (|TT| and |T2|) we must assume ^TZ — div uj — uP' 
and + div a; — a;^ to be constant on sections as well. This cannot be achieved unless 

div = (20) 

(on a sphere there are no vector fields of constant and non-vanishing divergence). We 
therefore have to assume 

-7^ - cu^ = const (21) 
2 ^ ' 

on sections. 

We may construct a metric go and a smooth vector field c^o satisfying (j2(Jj) and (j2H) in 
the following way: assume go to be axially symmetric 

go = d^2 + /(^)2 sin^^dy^l (22) 
/(O) = /(tt) = 1 
/'(O) = /'(tt) = 

and ijJq to be proportional to the axial Killing vector uj^ = A{9) d^. ujq is now automatically 
divergence-free and it vanishes at the two poles of the metric 9 = and 9 = n. 
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FIG. 1: Geometry of the horizon section 

In order to solve 1)211) assume that the curvature IZ attains its global minimum 'K-min at 
the poles, i.e. 71(0) > T^mm with equality at = and 6 = tt. Note that TZ is axially 
symmetric and therefore always has a local maximum, minimum or is constant near the 
poles. Moreover, all odd terms of its Taylor expansion at ^ = and 6 = tt vanish. Therefore 



A{9) 



7^.(0) TZmin (23) 



is well-defined and smooth everywhere. It also solves (PT|) . because \'Ti{6) — cuq ^ = |7^mm- 
The constructed horizon will have, in general, non-vanishing angular momentum with 
respect to (one can construct examples with strictly vanishing as well - see Appendix 
inj. It's interesting to note the geometry of sections is necessarily flattened on the poles 
(i.e. TZ is smaller) and rounded in the middle due to its rotation, just like in the case of a 
rotating drop of fluid, and in good agreement with intuition (see flg. 

We have to deal with the last equation now, pO|) . We need uj to rescale together with q in 
such a way that ()21|) is preserved. It's straightforward to see that we need the vector uj'^it) 
to be e~^^*^ujQ ^ (or uja^) = u;o a as a one-form). (piU) yields now an equation for GqAt": it 
must be proportional to oj 

GaAT'' =a{t)ujA. (24) 

Let's summarize the considerations now. We take the section metric to be g = e'^^^^^qo, 
the one-form ua = ujq a, where go and ujq are constant in time and constructed according 
to the prescription stated above. TZmin is the minimum of the curvature of q^, also constant 
in time. We introduce the short-hand notation G^^t^I^ = pi(t) and G^^r^k^ = P2{t). 
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(P)- (fT^ yield the following set of equations: 

F = ^UkO'^ (25) 
a = HkO'' (26) 
= -^nk-p,e'^ (27) 

9' = P2-'^n,e-'^-^n,9''. (28) 

All quantities except TZmin are functions of time. For a spacelike surface and TZmm > 
this system of ODE's can be reduced to the one in jl9|, which can be seen by substiting 
gF ^ ]ir^y2/3^ ^ -r2fe(p + ^), p2 = ni{p - ^), nkUi = -1/2 and eliminating 9^. 

We may obtain a Cauchy surface of an Einstein metric if we assume TZmin to be equal to 
(sections flat at the poles). We take Uk = pi = P2 = 9k = 0, F = const, ni(t) and nk{t) 
arbitrary. We can choose them in such a way that the horizon is spacelike. By comparing 
with ()B4|) and ()B3|) we conclude that in general the resulting initial data surface is neither 
time-symmetric, nor extremal. It is not a DH, as expansions of both null vector fields vanish. 
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APPENDIX A: GAUGE CONDITION FOR THE NULL ORTOGONAL VEC- 
TORS 

We assume that k'^ and are ortogonal to sections of H, which are preserved by the flow 
of n°. This is equivalent to 

^riJ'a ~~ ~l~ ^ ("^^) 

^n^a ~~ ^ ~t~ ("^^) 

for some functions A, B, C and D. Furthermore, the normalization conditions, differenti- 
ated, yield 

Vnka k'' = VJa = (A3) 
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VnKr = -VJak"". (A4) 

It follows that A = D = and B = —C. Without loosing generality and in order to simplify 
the equations we can put 5 = 0. This, substituted back to ()A2|) . leads directly to ((H). 

APPENDIX B: EXTRINSIC CURVATURE OF A SPACELIKE SURFACE 

For a spacelike surface (n/ > and < 0) we can normalize r and n 

r = (-2n,nfc)-i/V" = {-2nink)-^'^ n"^ . (Bl) 

The extrinsic curvature, i.e. the pullback of the convariant derivative of f, expressed in 
terms of introduced variables, reads: 

Kab = riiLAB + nkKAB (B2) 



ICajsn^'n'^ = \8nink\ ^^^i2„ln 

JCA/sn'^ = UJA + 7;Da In — 
2 rifc 



ni 



(B3) 
(B4) 



APPENDIX C: NON SPHERICALLY SYMMETRIC DH'S WITH VANISHING 
ANGULAR MOMENTUM 

We aim to construct a non-spherically symmetric dynamical horizon, of the family pre- 
sented in Section with strictly vanishing angular momentum. This can be achieved if we 
take f{9) such that 

1. 7^(7^/2 + 0) =n{n/2-e), 

2. 7^(0) > 7^™„ and 7^(0) = 7^(7^/2) = 7^(7^) = TZmin, 

3. all derivatives of TZ vanish at = 7r/2, 
and take uj = A{6) d^p, with 

A{e) = ± 



The positive sign is taken for < 7r/2 and negative for > 7r/2. A{6) is a smooth function 
despite the sign change, due to condition 3. is due to conditions 1 and 2 and ()C1|) . 
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Clearly, in this example the matter falling in does have angular momentum but its inflow 
is balanced to thanks to different directions of rotation above and below 9 — t:/2. 
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